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Abstract 

The mass of the glueballs is calculated in the B.S. equation frame- 
work. Under instantaneous approximation, the wave function of 
B.S. equations are obtained. The kernel is chosen as the sum of 
an one-gluon exchange potential, a contact interaction and a linear 
confining potential. The numerical results are in agreement with 
that of recent lattice calculation. 
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1. Introduction The existence of glueball states is a very important prediction of 
QCD, whose discovery would be a direct confirmation of the non-abelian (self-coupling) 
character of the gluonic degree of freedom in strong interactions. Therefore, these states 
are in focus of the interest of both theorists and experimentalists. Glueballs were 
suggested by Gell-Mann and Fritzsch[l|], Fritzsch and Minkowski 0. Since then a lot of 
investigation has been made. There are various theoretical approaches applied, which 
include the QCD sum rules P, |J, bag model [0, ||, nonrelativistic potential model [0] and 
lattice calculation]^, However, these approaches differ from each other markedly in 
their mass prediction, sometimes as large as 1 GeV. Therefore, it's necessary to make 
more investigations. 

On the other hand, in relativistic quantum field theory, the Bethe-Salpeter(B.S.) 
equation is an exact bound state equation. With some appropriate approximation it 
has been used successfully in the study of qq bound states[|10j- Therefore, it would be 



reasonable for us to try to study the bound states of the gluons in the framework of 
B.S. equation. In our investigation, the B.S kernel is chosen as the sum of an one- 
gluon exchange potential, a contact interaction and a linear confining potential. We 
calculate the mass of ++ ,2 ++ ,0 _+ and 2 _+ glueballs. The numerical results are in 
good agreement with that of recent lattice calculation. 

The paper is arranged as follows: in the next section, we construct the B.S. equation 
for the bound states of two constituent gluons, and then we discuss the integral kernel 
and the structure of the B.S. wave functions of glueballs investigated. In section three 
the numerical results and a discussion are presented. 

2. The B.S. Equation for Glueballs Let A®(xx) and A V ( X 2) be the gluon fields at 
points X\ and x%, \G) the bound state of two gluons with mass M and momentum P^. 
Then the B.S. wave function for a bound state is defined as 

X a ^(P,x u x 2 ) = (V\T{Al{ Xl )Al{x 2 ))\G), (1) 

where /x, v are Lorentz indices, and a, b color indices. The glueballs are color singlet 
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states, so we have 

X f u (P,x 1 ,x 2 )=5 ab x,u(P,x 1 ,x 2 ). (2) 
The translation invariance of system implies 

X^u(P, xi, x 2 ) = e lP ' x x^u(P, x), (3) 

where X = 1 — - and x = X\ — x 2 . We further define the B.S. wave function in 
momentum space as 

X,AP,q)= J ^ e -^X„(P,x), (4) 

where q is the relative momentum of the two constituent gluons. 

With a standard method we obtain the B.S. equation for a color singlet glueball 

state, 

X ,u(P, q) = A^( Pl )A^(p 2 ) J -^-G a ^(P, q, k) Xp a(P, k), (5) 



(2nY 

where the color indices have been suppressed. The tensor kernel G al3pa (P, q, k) in the 
above equation is defined as the sum of all two-particle irreducible graphs, and A^ a (pi) 
is the full propagator of the constituent gluons with momentum p iy and 

P = Pl+P2, (6) 

2q=p x -p 2 , (7) 

where P 2 =M 2 . 

Although this equation is formally exact, it is difficult to use in concrete problems. 
The reason is twofold: first we don't know how to calculate the kernel of B.S. equation, 
and second, even if we know the kernel we would be unable to solve the equation 
exactly. Therefore, to make the approximation is inevitable. 

In solving the B.S. equation of glueballs, we make "static" approximation as usually 
do in the case of qq mesons. In this approximation one neglects the dependence of the 
kernel G(P, q, k) on k and q . Then the dk integration on the right-hand side of 
equation (5) yields the equal time wave function 

c/V(P,k) = J dk oX ^(P,k). (8) 
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Next we perform dqo integration on both sides of equation (5), and at the left-hand 
side this again yields the three dimensional wave function </? Mi ,(P, q). On the right-hand 
side, the qo dependence is only contained in the variables pio and p 2 o- That is 

W(P, q) = J dq AMA^(p 2 ) J ^G a ^(P, q, k)^ CT (P, k). (9) 

To go on, an appropriate gauge has to be fixed in equation (9). It's convenient for us 
to choose Coulomb gauge. In Coulomb gauge the propagator is 



A^ik) — J2 ■ j- g^ u —(k^kv k^Sy k v s^) 



(10) 



k 2 + ie 
where 

s„ = (0,k). (11) 
We can write the propagator more clearly as 

A,(, ) = ^-( 5 ,-^), (12) 

Aoo(^) = (13) 
A 0l (k) = A l0 (k) = 0. (14) 

From (9) and (11), one finds that( for e M , the polarization vector of the constituent 
gluons, we set e = 0) 

Voi = = V^oo = 0. (15) 

Therefore, in Coulomb gauge one need only deal with the three dimensional wave 
functions (pij(P, q)(i = 1,2,3). It's appropriate to give an effective mass m to the 
constituent gluons and rewrite the propagator as 

A «« = FTTO'^ " & (16) 

For convenience, we adopt the center-of-mass system in the present paper. In this 
frame we have 

dq 2ni 



I (Pi 



m 2 + ie)(p 2 -m 2 + ie) E(4E 2 - M 2 )' 



where E = y/m 2 + q 2 . Substituting (16) and (17) into equation (9), we obtain 

£(M 2 -4£>y(P,q) = 

i(6w ~ q -f)(6 jr ~ qj f)j (03 G ^'( P ' q> k )^(^, k). (18) 
This equation is the starting point of our numerical investigation. 

The next question is how to construct the integral kernel. We assume that the 
kernel Gij k i(P, q, k) is the sum of two parts: the short distance part G\f k[ (P, q, k), and 
long distance part Gfj kl (P, q, k). As usual, we assume that the short distance part, 
G\f kl (P, q, k), is contributed by the three lowest order diagrams shown in figure 1. 
Besides equations (6) and (7), we also have 

P = P3+P4, (19) 

2k = p 3 -p 4 . (20) 
Calculating diagram a, b and c, G^} (P, q, k) can be expressed explicitly as 



G S*r( P ' k ) = 3 «( 47ra s) i^C tipT {p 1 ,p 3 )C uaT '{p2,PA) 

— {^9pv9p(J ~ 9pp9va ~ 9po9vp) j , 



„tO„t'0 „ri„r'j / / 

g 9 9 9 _ Mj \ 
p + /2 y 0i i \2 > 



(21) 

where /(= q — k) is the momentum exchanged between the two constituent gluons, and 
C^pApuPs) = (Pi - 1p?)p9pt + (pi + Ps)r^p + (pa - 2pi)p9pr, ( 22 ) 

C 1W r'(P2,P4) = (P2 - 2p4)v9(TT' + (P2 + Pa)t'9vo + {PA ~ ^P2)a9vr'- ( 2 3) 

The factor 3 on the right-hand side of equation (21) is the color factor which is | in 
the case of qq bound state. Diagram a and diagram b make the same contribution to 
physical states, so there is a factor 2 in equation (21). The strong coupling constant 
a s is chosen as a running one, 
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27 In^a+j^) 

QCD 



EZT' (24) 



where a is a parameter introduced to avoid the infrared divergence. As for the long 
distance part Gfj kl , we have no reliable knowledge about it, so we have to constructed 
it phenomenally. The experience with qq bound state shows that the long distance part 
and short distance part of kernel have different spin dependence. Generally speaking, 
the long distance part makes no contribution to spin-spin interaction. With the guid- 
ance of this experience, we only chose in equation (22-23) the terms containing tensor 
9iip9v<j as the spin dependence of the confining part because such terms have nothing 
to do with spin effect. Therefore we assume 

G$ p * = 2 *(Pi + Pa) ' (P2 + P4)9 w 9v*G(l), (25) 
where G(l) is spatial dependence of confining part of the kernel. We choose 

G(l) = 8 -f, (26) 

which corresponds to a linearly growing potential. The expression is very singular 
at the zero point of I, and regularization is necessary. The method is to subtract a 5 
function from the confining part of the kernel, that is to make the following replacement, 

8ttA 



G^lpa = 2i (Pl + Ps) ■ (P2 + P4)gpp9ua ^ 2 + 

- 5 3 (l) J d 3 k 1 2i(pi + p 3 ) • (P2 + PA)g w gv 
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(P + u 2 ) 2 ) 

where u is a small quantity. In actualcalculation we let u — > 0. In this way the infrared 
divergence is subtracted out. 

The next important question is to construct B.S. wave function for a given gluon 
bound state. For 0~ + glueball the most general decomposition of the four dimensional 
wave function is given by the following expression 

X ^(P,g) = / W PV. ( 27 ) 

For ++ glueball, 



Xpu(P, q) = /osw + h p p p v + f2 p p,qv + h p v% + hq^.q v , ( 2 §) 

5 



where / and fa are scalar functions of P 2 , q 2 and P ■ q. For other states the B.S. wave 
functions are more complicated and are not given here. Except for pseudoscalar state, 
the B.S. wave function has many independent components. This makes the numerical 
calculation difficult. However, when Coulomb gauge as well as the center-of-mass frame 
is adopted, the wave function becomes simple. 
From equation (18), we find 

q mj (P,q) = 0. (29) 

This equation gives a very strong restriction on the wave functions. The following are 
some three dimensional wave functions. 
For 0~ + state (pseudoscalar) 

(Pij(P, q) = f P {q)q k e ijk . (30) 

For ++ state(scalar) 



^■(P,q) = / S (q)(<%-^). (31) 



For 2 ++ state 



(Pij(P, q) = fti(q)(S ik - ^)r] k i(Sji - ^r) + /t2(q)(£*j - ^-)VkiQkqi- (32) 

For 2 _+ state 

<Pij(P, q) = f P i{(\)e k i m 'q mn q n (5 ik - < ^)(5 j i - -^) 

q z q z 

+ fp2(<i)(^k m nVmiqn ~ ^lmnf]mk1r) {$ik ~ ^t)(<^7 - ^f") ( 33 ) 
+ / P 3(q)?fc kVmnQmQn j 

where rjij is the polarization tensor of the glueball and f p , f s , f ti and f pi are scalar 
functions of |q|. 

Substituting equation (30) and the expression of kernel (21) and (25) into equation 
(18), we obtain an equation about 0~ + state(For simplicity, in the following formulae, 
|q| and |k| are expressed simply by q and k respectively.) 



E(NP - iE*)fJq) 



-8*f 



d 3 k 



k 

A(q, k)-cos6 + 8k 2 sin 2 6 
Q 



V s (q, k) + B(q, k^cosdViiq, fc)| f p (k), 

(34) 



where 6 is the angle between vectors q and k, and 

A( 4 , k) = W + f + V + 2, W - q2+ %-_%i s0 , (35) 
B(q,k) = M 2 + q 2 + k 2 + 2qkcos6, (36) 

47T 

K(<? ' *° = 3(q 2 + k 2 - 2qkcos9)ln(a + Z±*p^) ' (37) 

2A 

% ^ )= (^2^F' (38) 

The equations for other states are given in the Appendix. We can see that for 0~ + 
state and ++ state the equations are simple because each wave function only has one 
component. However, for 2 ++ state we have to solve two coupled integral equations. 
The appearance of the wave function of 2 h state is a little awful but the actual 
calculation shows that the three coupled equations about it are decoupled into two, so 
it is not too difficult to solve the equations. 

3. Results and discussion There are three parameters in the present model: a and 
Aqcd appearing in the running coupling constant and A, the string tension. Aq C d 
is chosen as 200 MeV which is commonly used in various models. The parameter a 
is chosen as 4.0. Such a choice implies that the running coupling constant tends to 
its largest value 1.0 when the exchanged momentum goes to zero. We can relate the 
parameter a to the threshold of the two massive gluon formation, and such a choice 
corresponds to a threshold of 400 MeV (2 Aqcd)- As for the string tension, A , it's 
different from the case of qq bound states where A is about 0.18 (GeV) 2 . As argued 
in ref. 0, in the most naive picture, the string tension between two gluons should be 
twice that between two quarks because each gluon acts like a qq pair. Therefore, we 
chose A = 0.36 (GeV) 2 . Such a value has also been used in glueball investigation [p]] 



and condensation calculation in QCD sum rules[12 



As most of the constituent models, in our investigation we give an effective mass to 
the constituent gluons despite the fact that gluons are massless in QCD Lagrangian. We 
can think that dynamical mass of gluons is generated through confinement interaction. 
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One measure of the effective gluon mass is the energy (« 2m) necessary to break the 
string joining two color-octet sources to materialize a gluon pair. Such a question has 
been studied in a lattice calculation with results m ~ 520 MeV. Non-perturbative 
continuum studies|n| also yield m = (500 ± 200) MeV. In our investigation we give 
the results when m = (0,400,600) MeV. 

With all the parameters determined, we solve the B.S. equation numerically. In 
actual calculation, the momentum integration must be cut at some value A, and when 
A large enough, the results are independent of it. The results are given in table 1. For 
comparison, the recent lattice results are also given in the table. 

We can see that our results are in good agreement with that of lattice calculation 
though there are some uncertainties due to the effective gluon mass. If we set the gluon 
mass to zero, the results favor that of the UKQCD group, while the results with a large 
gluon mass favor that of the GF11 group. 

In the early investigation, a light 0~ + glueball was preferred. For example, in MIT- 
Bag model§ Af(0 _+ )=0.4 GeV, in potential model0 M(0~ + ) = 1.4 GeV, and in 
QCD sum rules Q], M(0 _+ ) = 1.7 GeV. However, recent lattice calculation indicates 
that 0~ + is a very heavy particle |J, M(0~ + ) = 2.3 GeV. Our investigation gives the 
same result as that on lattice. 

In the nonrelativistic potential model, the form of the interaction between the two 
constituent gluons is similar to that of present paper, but their results are much sailer. 
There are two reasons for this fact.Firt, the B.S. equation is a relativistic equation. 
More relativistic effects still survive in our model, even many approximations have 
been made. Second, the string tension is an important parameter. In our model, the 
choice A = 0.36GeV 2 is reasonable. 

In fact, the effective mass of the gluon is still unclear. Someone has argued that the 
effective gluon mass is the result of confinement mechanism while the others think it 
incompatible with the principle of gauge invariance. It's interesting that our numerical 
results show that the glueball mass is quite insensitive to the effective mass of the 
gluons. We can even chose the gluon mass as zero without much change of the mass 
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of the glueball. 

It's known that nonrelativistic potential models are very successful in meson spec- 
trum calculation for heavy quark systems. One of the reasons of the success is that 
most mesons studied can be considered as nonrelativistic systems because the mass 
of mesons is quite near the sum of the mass of its constituents. However things are 
different for glueballs. Notice the fact, that the glueball mass is much larger than the 
sum of the effective gluon mass. This means that glueballs are relativistic systems. To 
deal with such systems, perhaps, the B.S. equation is a more suitable method. 

4. Acknowledgments. We wish to thank Prof. Yuan-Ben Dai and Prof. Xin-Heng 
Guo for interesting discussions. 

This work was supported by Chinese National Science Foundation. 
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Appendix 



The B.S. equation for ++ state 
E(M 2 -AE 2 )f s (q) = - 8 -/(^{ 



-A(q, fc)(l + cos 2 9) + A(q 2 + fc 2 )sm 2 ^ 



+ ±fl( ? , fc)(l + cos'Wfa, fc) - ^(3 - cos 2 ^)} 



V,(g,fc) 



(Al) 



The two coupled B.S. equation for 2 ++ state 



E(M 2 -AE 2 )f tl (q) = -Stt/HI 



-A(q,k)(l + cos 2 9y - k'sinW 



V s (q,k) 



~B(q, k)(l + cos 2 9) 2 V l {q, k) + ^(1 + cos 2 9) 2 } f tl (k) 
J d 3 k 



+8tt 



(2tt)* 



-A(q,k)k 2 sm 4 9 -2k A sin A 9 



V s (q,k) 



+^B(q, k)sin A eVi{q, k) + ^k 2 sin 4 9^ f t2 (k), 



E(M 2 -4E 2 ) f t2 (q) = 8n J 
k 2 



d 3 k 
(2^)3 



(A2) 



8q 



2 A(q, k)(5coc 4 9 - Qcos 2 9 + 1) 



2 {5coc 4 9 - 6cos 2 9 + 1) - 2{3cos 4 9 - 5cos 2 6> + 2) 



V a (g,fc) 



2g 2 

1 3 1 

+ ^B{q, k){5coc 4 9 - 6cos 2 9 + l)Vi(q, k) + ^jf (5cos 4 # - 18cos 2 9 + 5) j / t i(fc) 



-87T 



A" 1 



8q [ 



■A(q, k)(hcos 4 9 + 6cos 2 9 - 3) + 2£; 2 (3cos 2 6> - l)sm 2 # 



+ '—(5cos 2 9 - l)sin 2 9 
q 2 

3rv P 1 
+ ^^(5co S 4 ^ - lScos 2 ^ + 5) | f t2 (k). 



V a (q, k) + ^(g, k)(5cos 4 9 + 6cos 2 fl - 3)VJ(g, k) 



(A3) 
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The three coupled B.S. equation for 2 + state 

d 3 k 



E(M 2 — AE 2 ) (f pl ( q ) + f p2 ( q )) = -8n 



k 
q 

-8tt 



(2nf 



k 

-A(q, k)cos6 + 8k 2 sin 2 6 
Q 



V s (q,k) 



+ -cos6B(q, k)Vi(q, k) (cos 2 6 - -sin 2 9){f pl {k) + f p2 (k)) 



k 3 



d 3 k 



q 



A(q, k)cosOsin 2 e + 2k 4 sin 4 9 



V s (q,k) 



+— A (q, k)cos6sin 2 6Vi(q, k) j f p3 (k), 



E(M 2 - AE 2 ) (f pl (q) - 2f p2 (q)) = -8vr J 



+8k 2 sin 2 6 

„ f d 3 k 
-8tt / - — - 
J (2tt> 

k 3 



k 



d 3 k 



k 

—A(q, k)cos9 

q 



(AA) 



V s (q, k) + -cos9B(q, k^q, fc)| (f pl (k) + f p2 (k)) 
k 3 

— A(q, k)cos6sin 2 6 — Ak 4 sin 4 9 



V s (q,k) 



+ '— A (q, k)cos6sin 2 8Vi(q, k) j f p3 (k), 



E(M 2 -AE 2 ) f p3 (q) = -8rv J 
2k 4 



d 3 k 
~{2nf 



k 3 



(5cos 4 9 - Qcos 2 9 



k 3 

—A(q,k){5cos 3 9-3cos9) 
2q 



V s (q,k) 



(A5) 



(A6) 



+ ^cos 3 9 - 3cos9)B(q, k^q, k)j f p3 (k). 

where 9 is the angle between vector q and k, and the functions A(q, k), B(q, k), V s (q, k) 
and Vi(q, k) have been defined in equation(35-38). 
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Figuer Caption 



Fig. 1: The diagrams contributing to the short distance part of the B.S. kernel. 
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Table Caption 



Table 1: The mass(GeV) of the glueball states. 
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Table 1. The mass of the glueballs(GeV). 





M(0++) 


Af(0~ + ) 


M(2++) 


M(2-+) 


M(2++) 
M(0++) 


our results(m=0) 


1.50 


2.18 


2.17 


2.18 


1.45 


our results(m=0.4) 


1.60 


2.28 


2.23 


2.28 


1.39 


our results(m=0.6) 


1.75 


2.39 


2.32 


2.39 


1.33 


Lattice results [8] 


1.550± 0.050 




2.27± 0.100 




1.46 


Lattice results [9] 


1.740± 0.071 




2.359± 0.128 




1.36 



